arXiv:1505.04169vl [cond-mat.mes-hall] 15 May 2015 
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The Casimir interaction in a stack of equally spaced infinitely thin layers is investigated within the 
zero-frequency mode summation method. The response properties are considered to be described 
by a constant conductivity or by a Drude-Lorentz model with a finite set of oscillators consistent 
with the optical characteristics for graphite. It is found that the asymptotic distance dependence 
is affected significantly by the specific response. While the energy is ~ 1/d® for the constant 
conductivity model, the energy exhibits fractional dependence ~ 1/d®/® for the Drude-Lorentz 
description. The Casimir force on a plane is also strongly dependent upon the particular plane 
location in the stack. Furthermore, the calculated Casimir energy within the Drude-Lorentz model 
yields results in good agreement with measured cohesion energy in graphite. 

PACS numbers: 03.70.+k, 03.50.De 


I. INTRODUCTION 

Dispersive interactions originating from electromag¬ 
netic fluctuations play a prominent role for many ma¬ 
terials and devices. The Casimir force, first predicted 
by Casimir [1], is of importance especially for nanos- 
tructured or layered systems. One particular example 
is graphite composed of atomically thin graphene layers. 
The experimental realization of a single graphene [2] has 
attracted much attention due to its unusual properties 
and possible practical applications [3]. Contemporary 
experimental techniques enable the synthesis of not only 
a single graphene, but also a stack of a finite number of 
graphenes. Since graphene interaction is of fundamen¬ 
tal relevance, such progress in the laboratory motivates 
theoretical studies of how the electromagnetic coupling 
evolves as a function of the number of stacked layers. 

The graphene/graphene and graphene/substrate 
Casimir interactions have been studied by several 
authors recently. Earlier reports have utilized the 
hydrodynamic model for thin conducting shells [4], 
however these results were found to be not suitable since 
the hydrodynamic model does not take into account 
the Dirac-like nature of the carriers [5]. Other studies 
have described the graphene response properties via the 
longitudinal polarization function, optical conductivity, 
and the polarization operator in a 2 -|- 1 model [6-8]. 
These approaches have solidified our understanding that 
the Dirac-like carriers coupled with other factors, such 
as temperature and doping, result in highly unusual 
behavior of the Casimir force in terms of distance 
dependence and other asymptotic dependences [9]. 

The optical conductivity of graphene a is of particular 
interest for the Casimir interaction. It has been shown 
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[10] that a takes a constant universal value {agr = e^/Ah) 
over a relatively large frequency range, hoj < 3 eV. This 
frequency range is commensurate with distance separa¬ 
tions in the /im and sub-/iTO scale, meaning that agr 
is suitable to determine the asymptotic behavior of the 
graphene Casimir interaction [7]. Considering other con¬ 
figurations with an arbitrary value for the constant con¬ 
ductivity has also been of interest, including calculations 
for the Casimir energy for planar and spherical symme¬ 
tries [11]. The energy dependence on the separation be¬ 
tween two planes is £ ~ A/dP, while for a sphere with 
radius i? it is £ ~ 1/R. One arrives at these results by 
simple dimensional considerations. In both cases, taking 
the conductivity a to be the graphene universal value, 
agr, leads to a Casimir interaction that does not depend 
on the Planck constant and the velocity of light in agree¬ 
ment with earlier results for graphene planes [7]. 

The purpose of this paper is to investigate the 
vdW/Casimir interaction in a stack of Af infinitely thin 
parallel layers, as each layer is described by its con¬ 
ductivity a. Different models are taken for a including 
a = const, which can be appropriate for graphene sheet 
and a = oo, which is suitable for perfect metals. Results 
for a Drude-Lorentz (DL) model are also derived with 
the application for graphite parametrization [12]. 

We utilize the mode summation method with the zeta- 
regularization approach for the calculations of the inter¬ 
action energy and force. The main ingredients of this 
theory with relevant applications are summarized in [13]. 
The energy calculated via the mode summation is a reg¬ 
ularized quantity in the following form 
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where 4' defines the energy spectrum from the relation 
4'(a;) = 0. The parameter A with a wavenumber dimen- 
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sion is introduced to preserve the energy dimension of 
£’(s). Note that £(s) must be renormalized by subtract¬ 
ing corresponding to the energy for infinitely sepa¬ 
rated objects (d —>■ oo) 
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The paper is organized as follows. In Sec. II we derive 
the general expression for Casimir energy for a stack of 
Af planes with constant conductivity and analyze them 
in detail. Section III is devoted to the derivation of the 
Casimir force acting on a plane in the stack with constant 
conductivity. In Sec. IV we consider the DL model with 
parametrization for graphite and analyze expressions ob¬ 
tained for the energy. We show that the binding energy 
per unit atom is close to the one obtained experimen¬ 
tally. Section V contains our conclusions and discussion. 
In Appendix A we make derivation of the expression for 
the energy of stack of Af planes. The calculation of the 
force is made in Appendix B. 

Throughout the paper we use units in which h = c = 1. 
Where necessary the constants are restored for better 
understanding of the final results. 


II. CASIMIR ENERGY IN A STACK OF 
PLANES 


interaction energy stored in the stack can be repre¬ 
sented as 

1 pOO pi 

J y'^dy J dx In Fj^ (^r]x, , (4) 

where rj = jc is a dimensionless parameter and the 
integrand function is 
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where / = f{t, z) = (cosh z Ft sinh zY — I -I- coshz -I- 
t sinh z. Details of the derivation of the interaction energy 
can be found in Appendix A. It is noted that the cases of 
A/" = 3,4, consistent with Eq. (3), have been examined 
in Ref. [14]. 

To understand the main factors affecting the Casimir 
energy, we define the average energy per plane for the 
stack of planes of N planes 


The system under consideration consists of N equally 
spaced parallel and infinitely thin layers, as shown in Fig. 

1 . 



FIG. 1. A schematic representation of N parallel infinitely 
thin layers separated in the 2 -direction and extending in the 
X — y plane. The planes are considered to be identical and 
equally spaced according to dm = md. 

We first calculate the Casimir energy for this system 
by assuming that all planes have a constant conductivity. 
For this purpose, the appropriate boundary conditions 
for the electric and magnetic fields are imposed, which 
leads to a separation of the transverse electric (TE) and 
transverse magnetic (TM) modes [11]. As a result, the 


^(Y) ^ ^ ^ Q'fu W) + Q^TE\y) /.A 

Af d^ ^ ^ 

where Q^'^\'r]) = /Af. 

We further consider the energy ^ for infinite num¬ 
ber of planes, which can be expressed by taking the 
A/" —> 00 limit in Eqs. (3) - ( 6 ). As a result, the function 
Fjy{t,z) is replaced by $(t, z) = limA^^oo z) = 

e"V(C 2 )/(l +0 giving Qt'e.tm- 

A. Calculation of the Casimir energy 

It is instructive to study the case of infinite number of 
perfectly conducting layers firts. We see that for 77 00 , 

one finds Fy (1 ~ from Eq. (5). Therefore, 

the energy stored in the stack is 

£:W) = (Af-l)gcas, (7) 

2 

where Scas = ”75033 is the energy between two perfect 
metallic substrates separated by a distance d. Using Eq. 
( 6 ) for the average energy, it is easy to see that for infinite 
number of planes, we recover the standard result for the 
Casimir energy between perfect metals: 

lim = £cas- ( 8 ) 
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In addition to ry —>■ oo, we also calculate the energy 
for a small constant conductivity in the limit of ry —> 0. 
In this case, the main contribution comes from the TM 
modes, and the energy per unit plane is of the form 

= ScasabM, (9) 

where 

(;■!)■ 

These expressions give the means to evaluate the 
Casimir energy for different number of layers M. For 
example, for two layers (M = 2) and infinite number of 
layers {M = oo), in the 7y —>■ oo we find the analytic 
relation 

62 = 1.4864; &00 = ^. (II) 

Numerical calculations for intermediate cases are pro¬ 
vided in Table I. It is clear that the bj^ coefficient plays 
the role of a scaling constant to the energy. Therefore, 
the interaction has a universal behavior, in which the en¬ 
ergy per unit plane is simply proportional to £cas and 
the small conductivity. 


AT 

bM 

N 

bjsf 

2 

1.48641 

5 

2.46766 

3 

2.02535 

6 

2.57937 

4 

2.30086 

OO 

TT 


TABLE I. The constants jEcasO in the case of 

small conductivity cr —> 0. 

Numerical simulations help us understand how the en¬ 
ergy behaves for intermediate values of ry. In Fig. 2 
results are shown for the energy per unit plane normal¬ 
ized by Ecas for different number of the planes. We ob- 





FIG. 2. The normalized average energy of M planes 
jJsfScaa as a function of the dimensionless conductivity 
jy = 2-eg jc. 


serve that increasing M correlates with larger energy. For 
infinite number of planes the energy per unit plane ap¬ 
proaches the Casimir energy for two ideal metals. 


III. THE FORCE 


To gain a better perspective of the Casimir interaction, 
the force on a particular plane in the stack M is also con¬ 
sidered. This is a convenient way to explicitly understand 
how the interaction changes depending on the location of 
the plane. The Casimir force acting on plane m in the 
stack M (Fig. I) can be written as (see Appendix B) 










( 12 ) 


where the TE and TM contributions read 

^ (L I). 

UtT'^\v) = J y^dy J dxGm.Af (yx, . 

These are expressed in terms of a single function 

_ J2(Af-l)) _ (1 _ + t) ’ 

( 13 ) 

where /(t, z) was defined in the previous section. Ex¬ 
amining the general properties of Eq. (13), one finds 
the following. The function G is odd with respect to its 
indexes, thus Gm,M = —Gjy-i-rn,M- Therefore, for two 
symmetric planes, for which m Af — 1—m, the force has 
the same magnitude, but opposite sign signaling attrac¬ 
tion. Also, the force on the central plane m = {M — I)/2 
in a stack of odd number of plane is zero. 

Perfect Metals ry —>■ 00 . We first consider the case of 
perfect metallic planes. The force on a plane has different 
characteristic depending on the particular m. Eor the 
first plane at m = 0, it is obtained that lim,;_,.oo Go, tv/ = 
— , which leads to equal TE and TM contributions 
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Therefore, the force on the very first (or last) plane and 
the rest of the stack corresponds to the Casimir force 
between two perfect metals. 

For the second plane m = I in the stack, the TM 
contribution is: 
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For the TE modes, however, the domain of integration 
over X needs to be divided into two parts to avoid the 
occurring divergences. Thus it is found 


lim = lim f y'^dy [ dxGi^f^ (yx, 
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Further examination shows that the first term in the 
above equation falls down as r]~^ and second term - as 
77 “^. Comparing with lim^_>oo ~ if is con¬ 

cluded that the dominant term is the TE term of the 
order r]~^. Nevertheless, the force on the second plane 
will be asymptotically zero: = 0 . 

For any other plane Af/2 > m > 2, the TM contribu¬ 
tion is evaluated by changing variables y = ln(l -|- sx/if) 
in Eq. (12), after which the limit 77 00 is executed. 

One obtains 

g3^(^2m _ 

_ 2^ 

( 17 ) 

where w = 1 -I- s -I- 1 / 5 ^ 5 ”+^. The results for the TE 
modes are the same as for the m = 1 plane (Eq. (16)). 
We see that for any plane with m > 2 has even 

stronger dependence on the conductivity as compared to 
the second one. Nevertheless, dominant term in is 

of the order r]~^ (Eq. (17)) and it ensures that the force 
on all planes except the first (or last) one is asymptoti¬ 
cally zero. This is a curious result, which says that the 
dominant interaction present in a stack of perfect metal¬ 
lic planes is between the outer one and the rest of the 
stack with the Casimir force being the typical one for 
two perfect metallic substrates. 

Small conductivity y —> 0. The case for vanishingly 
small conductivity is considered next. For this situation, 
the force on the m-th plane is determined primarily by 
the TM mode according to: 

pOO poo / -I \ 
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Evidently the force on any plane in the stack will have 
the same distance dependence and proportionality to a. 
Its strength is determined by the numerical value of the 
above integral for a given m and M. To quantify how the 
particular location of the plane affects its Casimir force, 
we define the following parameter 

— y: ) (79) 

Cas 

where iFcas = A/2A0d‘^ is the Casimir force for two per¬ 
fect metals. For Af = 00 and large m, we obtain the 
following asymptotic expression 


,00 


2880(ln8-2) 
A{2m + 1)4 


( 20 ) 


showing that the Casimir force falls down quickly as the 
plane moves away from the end of the stack. This ~ m~^ 
behavior can be compared with the numerical calcula¬ 
tions for Sm»i,oo using Eq. (18). The first 10 coefficients 
are shown in Table II for the case of the infinite stack. 
We see that the asymptotic expression Sm,oo differs by 
mere 7% for m = 2 from the computational result, while 
for the seventh plane the error is smaller than 1 %. 


m 

^m,oo 

^Tn,oo 

m 

^Tn,oo 

^m,oo 

0 

3.05 

14.51 

5 

4.9- 

10-4 

1.5- 

10-3 

1 

7.6- 10"^ 

2.4- 10"^ 

6 

2.5- 

10-4 

7.7- 

10-4 

2 

1.1 • 10"^ 

3.3- 10"^ 

7 

1.4- 

10-4 

4.3- 

10-4 

3 

2.9 • 10"^ 

9.1 ■ 10"^ 

8 

8.7- 

10-® 

2.5- 

10-4 

4 

1.1 • 10"^ 

3.3- 10-3 

9 

5.6- 

10-® 

1.5- 

10-4 


TABLE II. The constants Sm,oo = /J-casix in the 

case of small constant conductivity ct —>■ 0 and Sm 00 = 
A'^'^\nDL)/Fcasa,r for the Drude-Lorentz model. 


Numerical Calculations For intermediate values of 77 , 
one must examine Eqs. (12) via numerical techniques. In 
Fig. 3, we show how the normalized to the perfect metal 
value force evolves as a function of the conductivity in 
the case of Af = 00 . 


Jim- 



FIG. 3. The Casimir force normalized by iFcas = 

2 

as a function of the dimensionless conductivity y = 2TTa. 
The inserts show the normalized Casimir force vs. 77 for the 
m = 1 and m = 2 planes. 

It is observed that the force starts out as ~ 77 and then 
asymptotically approaches the value for perfect metals. 
The force on any other plane {jn = 1,2 only shown), 
however, exhibits an extreme behavior with a sharp max¬ 
imum in the small 77 range. Also, for any ttt > 1 

is orders of magnitude smaller than showing that 

the overall Casimir interaction is determined primarily 
by the very first plane and the rest of the stack for any 
value of the conductivity. 


IV. DRUDE-LORENTZ MODEL OF 
CONDUCTIVITY 

The Casimir interaction depends strongly on the re¬ 
sponse properties of the materials. In addition to in¬ 
finitely thin planes with a = const, we also calculate the 
Casimir interaction within a Drude-Lorentz (DL) model 
for the conductivity of each plane. The DL model leads 
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to an interesting asymptotic distance dependence show¬ 
ing once again the richness of this ubiquitous interaction. 

As a model, we utilize the optical response of graphite. 
First principles calculations indicate that the optical con¬ 
ductivity as per graphene sheet is very close to the one 
for an isolated graphene over a wide range of frequen¬ 
cies [15]. The results for graphene have been mapped to 
a Drude-Lorentz model consisting of a Drude term and 
seven Lorentz oscillators according to [12] 


o-dl(w) 


/o^p ^ ^ 

7o — iw 

j — \ J 


( 21 ) 


In this case, the Casimir interaction is calculated using 
Eq. (2), where the spectral function is a function of 
ctdl, tT(w)). It is realized that 

9,,; In'E'= ;^ (9,^'i;'-I-, (22) 

has the same simple pole structure as The renor¬ 
malized expression d^j In , however, has no poles. 

Therefore, the Casimir energy can be found using Eq. 
(2) together with the formalism in Sec. II by setting a 
to the DL conductivity. 

The graphene DL conductivity is obtained from Eq. 
(21) for the 3D graphite by multiplying it with 27ra/c (a 
is a characteristic distance typically taken as the inter- 
planar distance of graphite). The expression is given on 
the imaginary axis uj = iXc (k = iX) as follows: 


Vdl{X) 


Volo Xrik% 

70 + A ^ A2 -t -b Ayfe ■ 


(23) 


Here, 7 fc is the relaxation time and ujk is the characteristic 
frequency for the fc-th term. Also, jk = Jk/c, Xk = 
Wfc/c, and = 2TTafhU!p/cjk- In Table III we show the 
parameters using the calculated values for graphite [12]. 


rjk 


7fe 

eU 

OJk 

eV 

% 

1 

nm 

Afc 

1 

nm 

Vo 

0.01712 

70 

6.365 

- 

- 

70 

0.0322 

- 

- 

Vi 

0.13855 

7i 

4.102 

LJ-l 

0.275 

71 

0.0207 

Ai 

0.0014 

V2 

0.05949 

72 

7.328 

U}2 

3.508 

72 

0.0371 

A 2 

0.0177 

73 

0.37991 

73 

1.414 

0 J 3 

4.451 

73 

0.0072 

A 3 

0.0225 

74 

0.08462 

74 

0.46 

6J4 

13.591 

74 

0.0023 

A 4 

0.0688 

75 

1.09548 

75 

1.862 

tos 

14.226 

75 

0.0094 

As 

0.0721 

76 

0.30039 

76 

11.922 

LOG 

15.55 

76 

0.0604 

Ae 

0.0788 

77 

0.03983 

77 

39.091 

LOj 

32.011 

77 

0.1981 

A? 

0.1622 


TABLE III. Parameters for the Drude-Lorentz model of a 
graphene sheet in 3D graphite 


We note that the optical response in the infrared 
regime for 3D graphite and an isolated graphene is differ¬ 
ent. While cr for graphite exhibits a Drude-like behavior, 


the graphene optical conductivity is constant. This dif¬ 
ference is attributed to the different electronic structure 
characteristics for the two systems in this range [15]. To 
ensure that the a = const feature is captured, the single 
graphene conductivity rjjjL is obtained by using a char¬ 
acteristic distance a = 0.224 nm. In addition, we require 
that ?)dl(0 ) coincides with rjgr at zero frequency as: 

VDLiX) = 7dl(A) —. (24) 

Vo 

The Casimir interaction energy in a stack of planes 
characterized with a DL response can now be calculated 
using Eqs. (3) - (6), but with conductivity that is de¬ 
pended upon distance - 7 (If)- Examining the large 
separation limit d —>■ oo {qj = 2djj oo) shows that 
tjdl Vo {vdl Vgr)- Therefore, the Casimir energy 
has the same behavior as the one described in Sec. 11. 

For small separations d —>■ 0 {qj —>■ 0), where the 
planes approach the formation of graphite, we find rj « 
if Vjlj 0- The main contribution to the energy 

comes from the TM mode. Using Eq. (3), we obtain 




r 

327r2 


c(Af) 


(vdl) 


Rn' 

d5/2’ 




(25) 


where r = hvjljj ^ parameter collecting the 

characteristic frequencies and oscillator strengths from 
the DL conductivity. The d“®/^ asymptotic distance de¬ 
pendence is a consequence of the frequency dependence 
of the response properties. We also realize, that all os¬ 
cillators in the DL model participate in the interaction 
via the r, which reflects the broad band nature of the 
Casimir phenomenon. 

Therefore, as the planes approach the large and small 
separation limits, the energy exhibits a transition pro¬ 
cess in terms of the asymptotic distance dependence. To 
understand this better, we consider a system from Fig. 
1 with N = 2. Starting out at d —)■ oo, where the planes 
behave as graphene sheets with ry = const, the energy 
£ ~ d~^. As the planes become very close and approach 
the graphite interplane separation a — 0.334 nm, the con¬ 
ductivity is described by the DL model. Thus the energy 
exhibits £ ^ d“^/^ asymptotic behavior. This transi¬ 
tion is shown in Fig. 4 representing plots of the ratio of 
the Casimir energy for two planes with DL conductivity 
rjDL and the Casimir energy for two graphene planes with 
constant conductivity rjgr- The maximum of the energy 
ratio at d « 14 nm indicates that as d increases, the DL 
nature of the response becomes less relevant. 

Note that a fractional d-dependence of the Casimir in¬ 
teraction was obtained in other systems as well. For ex¬ 
ample, the energy between two infinitely thin metals has 
the d“®^^ characteristic behavior as found by ab initio 
simulations [16]. Thin metallic films, described by a 
Drude model dielectric response, are also characterized 
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fi<"’('7DL) 


S‘"'('7gr) 



d(nm) 


FIG. 4. The ratio of the Casimir energy for two planes with 
dimensionless DL conductivity tidl and the Casimir energy 
for two graphene planes with constant dimensionless conduc¬ 
tivity i]gr. In the framework of the DL model, the energy falls 
down as d~^ for large separation and as for small sepa¬ 

ration. The graphene interaction energy has dependence 
for any distances. 


by a fractional d-dependence [17]. Based on these re¬ 
sults and our investigation, we may conclude that for 
nanostructured systems described by a DL response, the 
Casimir interaction is expected to have an unusual, frac¬ 
tional asymptotic distance dependence. 

We further consider the average energy per unit plane 
in an inhnite stack of planes using Eq. (3). Our re¬ 
sults show that the separation a plays an important 
quantitative role for £. While a = 0.334 nm for the 
DL model is used for the graphene planes in graphite, 
a = 0.224 nm is used to ensure that when graphene 
planes are separated far apart, the graphene conduc¬ 
tivity approaches its constant value r]gr. We obtain 
that when the plane-plane separation is d = 0.334 nm, 
^{iIdl) = 454.2^, £{rjDL) = 357.1^ for the two 
cases, respectively. Describing each plane with a constant 
conductivity yields £ = 65.2;^ for the same separation. 
Therefore, the average energy obtained via the DL model 
is « 5 — 7 times greater than the one for rjgj.. 

Estimates for the graphene binding energy in graphite 
£ib = £lap can also be given. Here p = 2.83g/cm^ 
is the graphite mass density and a = 0.334 nm 
is the graphite interplane separation. We find that 
£ibiVDL) = 74 meV/atom, £ib{^DL) = 58 meV/atom 
and for constant graphene conductivity model £ib('r]gr) = 
11 meV/atom. First principles studies report cohe¬ 
sion energies 24 -L 26 meV/atom [18], 24 meV/atom 
]19]. Experimental data shows that the cohesion is 
35±10,15 meV/atom [20], 61 ±5 meV/atom [21]. There¬ 
fore, our results for graphite sheets described via a con¬ 
stant conductivity underestimate the binding energy, 
while the DL response ones are fairly in line with the 
experimental data. 

To understand how the interaction evolves as a func¬ 
tion of the particular DL plane m in the stack (Fig. 1), 
the Casimir force is considered using Eqs. (12), (13). 


The Sm,N' from Eq. (19) with a = agr not only depends 
on m but also on the interplanar distance d, which is 
unlike the constant conductivity case. Furthermore, the 
force on the m-th plane is a function of the combination 
(2m ± 1 )? 7 , which can be used to derive asymptotic ex¬ 
pressions in terms of the m-th dependence. It is found 
that for {2m+l)fjmax ^ 1 witht m > 3 and fjmax = 0.177 
corresponding to the maximum of the DL conductivity, 
the following expressions hold: 


~TM 


^TE 


37800 

T^^V?naxVgr{‘2m + 1)^ 
1200 

r]gr{2m + 1)6 



Therefore, the large m behavior is determined primarily 
by the TE mode contribution according to Eqs. (26). 
To compare with the results for small m, we give the 
first 10 values of Sm,oo in Table IT It is evident that 
the magnitude of decreases rapidly beyond the first 
few planes in the stack, similar to the case of constant 
conductivity. 


V. CONCLUSION 

In summary, we have applied the zero-frequency mode 
summation method to calculate the Casimir interaction 
in a system composed of infinitely thin equally spaced 
layers. The obtained boundary conditions lead to sep¬ 
arating the TE and TM mode contributions, which is 
typical for planar configurations. The response proper¬ 
ties are considered within a constant conductivity and 
Drude-Lorentz models. 

The energy stored in the stack of planes characterized 

with constant conductivity is found to be ^ 1/d^ 

with a magnitude dependent on a and the number of 
planes. While for perfect conductors, the usual Casimir 
energy for two infinitely conducting metals is recovered, 
^ a for small conductivity planes. Interestingly, the 
Casimir force has a very strong dependence on the loca¬ 
tion of the plane in the stack. The attraction is strongest 
between the first few planes and the rest of the stack and 
it decreases by orders of magnitude for larger m. 

The DL model changes dramatically the asymptotic 
distance dependence of the interaction especially at layer 
separation d < 5nm. The energy has a fractional depen¬ 
dence according to £^'^^ ~ l/d®/^. The characteristics 
of all oscillators, quantified in the parameter r, enter the 
energy expressions indicating that all oscillators are im¬ 
portant. Also, using the specific values for the graphite 
DL conductivity, we are able to obtain binding energies 
consistent with reported experimental data. 

This study broadens the application of the zero-point 
summation method beyond the typical Casimir interac¬ 
tions between two objects of planar symmetry. It shows 
that this particular approach can be used not only in 
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systems with constant conductivity, but also in systems 
with frequency dependent response. Our results present 
further evidence of the complexity of this phenomenon in 
terms of the nontrivial effects and importance of geome¬ 
try and response of the interaction objects. 
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Appendix A: Interaction energy stored in the stack 
of J\f planes 


There are A/” — 1 interplanar regions for the stack of 
Af planes, and the vacuum modes in each region are de¬ 
fined by two constants for a total of 2{N' — 1 ) constants. 
The spectrum outside of the stack is defined by a sin¬ 
gle constant in each semi-space. Therefore, there are 2JV 
constants to be defined. At the same time, we have 2JV 
equations for the modes (see Ref. [11]), which can be 
found from the appropriate boundary conditions on a 
plane at a position z = dj for the TM mode 


K]z=d, = 0 ^ 




Aziia , . 


and the TE mode 


(Al) 


[^z\z=dj ~ 

[*z]z=d, =47rfcra;/i^|^=d,, (A2) 

with the electromagnetic fields being 




Here e{z),h{z) are the fields amplitudes, which depend 
on the z coordinate only and [f{z)]a = f{a—0) — f{a+0). 
Each type of mode is considered separately below. 

TM mode. The main determinant of the system of 
equations Eq. (Al) is expressed in terms of the following 
matrix 


/ _eP. e-w \ ^ ^ ( -1 0 \ 

y(-6-l)eP^ (6-l)e-P^y’ ^ \^-6-10y’ 






-e-w-i 

e-Pj-i 


1^0 


(A4) 


where pj = djK, k = + ky + A^, and w = iX. Also, 

b = 2rjK/\ and 77 = 27rcr (= 27r(T/c in dimensional vari¬ 
ables) . The determinant has the following Af x Af block 


matrix: 


A = 


0 

0 


Bi 0 
Al 82 
0 A 2 


0 0 0 
Z 2 0 0 


0 0 

0 0 

0 0 

Aj\f-2 

0 Aj\f-i 


(A5) 


The determinant is then transformed in a diagonal 
block form 


Z 3 0 ... 0 0 

0 Al ^2 0 0 

0 0 A 2 ... 0 0 

A= : : : . . , 

0 0 0 ... Aj ^-2 Bj ^-1 

0 0 0... 0 Aa^_i 

where 

n—1 

Z,=Zi-{-irl[{B,A-^)Z2. 

1=1 

This enables writing A in the following way: 

A = 2^-1 det [Z 3 ], 


(A 6 ) 


(A7) 


(AS) 


where we have taken into account that det Aj = 2 . This 
factor is omitted in what follows as it does not change 
the energy spectrum. 

Taking into account that the planes are equidistant, 
the determinant can be cast a slightly different form: 


A = det 


Al 




where p = dn and 


(A9) 


G = B,Ay 


-1 


C — cosh p — b sinh p sinh p 
5 cosh p -I- sinh p — cosh p 


(AlO) 


The Af — 1 power of the matrix C is executed by a Jordan 
form representation C = TJT~^ with 


C coshp+/_ coshp+/+ \ 

sinh p+6 cosh p sinh p+fe cosh p | 

1 1 J 




where f± = ±y (coshp J-| sinhp)^ — 1 — (coshp -|- 
I sinhp). Therefore, we arrive at 


A = 


g-p(N-l) 



1 _ pM-2 

l-P 


l + |l-/^A \ 

/ I-P )' 


where / = |/_| = y (coshp J- | sinhp)^ — 1 -|- coshp -|- 
I sinhp. 

















To renormalize the spectrum, the limit d —> oo (p —>■ 
oo) is calculated resulting in 


Aq = lim A = — 

p—VOO 



(All) 


The renormalized determinant reads 


/b \ A 


X 



1 _ / 2 (A^- 1 ) 

I- P 


1 + 1 l-/2Ar \ 

/ I- P )' 


(A12) 


To obtain the energy (see Eq. (2)) the following inte¬ 
gration is needed 


Appendix B: Force on the m-plane 

The Casimir force upon a particular plane m can be 
calculated using Eqs. (A 8 ), (A14) and (A16) by taking an 
appropriate derivative with respect to dm- For this pur¬ 
pose, we consider the following ratio lime_j.o(ln A((m -I- 
e)d) — In A{md))/e. The matrix can be regrouped by 
realizing that dm appears in Am and Bm+i only. 


Z3 = - (- 1 )-^ 

[i=i 

X {iBmAm^){l3m+lAm\l)} 

x| n (BlApz,. 

I j=m-|-l I 


(Bl) 


f /■°° dkj;dky 

JL^ 


dX In Fj^. 


(A13) 


Changing the integrand variables to spherical coordinates 
kx = Ksindcosp,/cy = Ksindsinp, A = kcosB (k = 

^kx -I- + A^), y = 2p = 2dK,X = kx,x = cos 9 (note 

that b/2 = rj/x) the TM contribution to the energy be¬ 
comes 


(,(At) _ QtM id) 

^TM — ^3 > 


(AM) 


where 

Q™(d) = j^^ (x’ I) ■ 


and the dimensionless function F^(t,z) is given in Eq. 
(5). 

TE mode. The spectrum calculation for the TE mode 
follows the same path with the substitution b = 2rjK/\ —>■ 
2r}X/K. Thus the TE energy is found as 


where 


(.(At) _ QT^^(b) 
TE - dS ■- 


(A16) 


Qt'e^(?7) = 


327r2 


) pi 

y'^dy / dxlnEy {r]X, . (A17) 


Taking dm = (™ + e)d, one can expand in first order of 
the small parameter e 

iBmAp)iBm+lAp^^) = fj +ebp . (B2) 

Thus we can write the following relation 

A((to + e)d) = det(M + ebpP) = A det(l -|- ebpM~^P) 
= A(l + e 6 KTr(M-ip)). (B3) 


where 

At-l 

M = - (-1)-^ n {B,A-^)Z2 

f=i 

= Zi - 

1 —1 — m ry 

b -p 

and the matrix C is calulated in (AlO). Thus we find 

lim ^ bKTr{PM-^) 

€->■0 e 

1 K62/e-P (/2™ - /2(At-l-m)^ 

= 2 e-P/(l - /2(A-1)) _ (1 _ /2At) (1 + I) ■ 

In the limit of d —> oo (p —)■ oo) the above expression 
tends to zero meaning that there is no need for renor¬ 
malization. Changing the integrand variables to spheri¬ 
cal coordinates as previously done in Appendix A, Eqs. 
(12), (13) are obtained. 


P = -{-ipc 


[1] H. B. G. Casimir, Proc. K. Ned. Akad. Wet. 51, 793 
(1948); H. B. G. Casimir and D. Polder, Phys. Rev. 73, 
360 (1948). 

[2] A. K. Geim and K. S. Novoselov, Nat. Mater. 6, 183 


(2007). 

[3] M. Katsnelson, Graphene: Carbon in Two Dimensions 
(Cambridge University Press, 2012) p. 351. 

[4] G. Barton, J. Phys. A: Math. Gen. 37, 1011 (2004); J. 














9 


Phys. A 38, 2997 (2005); M. Bordag, J. Phys. A 39, 6173 
(2006); Phys. Rev. D 76, 065011 (2007); M. Bordag and 
N. R. Khusnutdinov, Phys. Rev. D 77, 085026 (2008). 

[5] G. L. Klimchitskaya and V. M. Mostepanenko, Phys. 
Rev. B 91, 045412 (2015). 

[6] B. E. Sernelius, Phys. Rev. B 85, 195427 (2012). 

[7] D. Drosdoff and L. M. Woods, Phys. Rev. B 82, 155459 

( 2010 ). 

[8] M. Bordag, I. V. Fialkovsky, D. M. Gilman, and D. V. 
Vassilevich, Phys. Rev. B 80, 245406 (2009); I. V. Fi¬ 
alkovsky, V. N. Marachevsky, and D. V. Vassilevich, 
Phys. Rev. B 84, 035446 (2011); M. Bordag, G. L. Klim¬ 
chitskaya, and V. M. Mostepanenko, Phys. Rev. B 86 , 
165429 (2012); M. Bordag, G. L. Klimchitskaya, V. M. 
Mostepanenko, and V. M. Petrov, Phys. Rev. D 91, 
045037 (2015). 

[9] G. L. Klimchitskaya, U. Mohideen, and V. M. Mostepa¬ 
nenko, Phys. Rev. B 89, 115419 (2014); G. L. Klimchit¬ 
skaya and V. M. Mostepanenko, Phys. Rev. A 89, 052512 
(2014); G. L. Klimchitskaya, V. M. Mostepanenko, and 
B. E. Sernelius, Phys. Rev. B 89, 125407 (2014). 

[10] V. P. Gusynin, S. G. Sharapov, and J. P. Garbotte, 
J. Phys.: Gondens. Matter 19, 026222 (2007); L. A. 
Falkovsky and A. A. Varlamov, Fur. Phys J. B 56, 281 
(2007); R. R. Nair, P. Blake, A. N. Grigorenko, K. S. 
Novoselov, T. J. Booth, T. Stauber, N. M. R. Peres, and 


A. K. Geim, Science 320, 1308 (2008). 

[11] N. Khusnutdinov, D. Drosdoff, and L. M. Woods, Phys. 
Rev. D 89, 085033 (2014). 

[12] A. Djurisic and E. Li, J. Appl. Phys. 85, 7404 (1999). 

[13] M. Bordag, G. Klimchitskaya, U. Mohideen, and 
V. Mostepanenko, Advances in the Casimir Effect (Ox¬ 
ford University Press, Oxford, 2009) p. 749. 

[14] N. Khusnutdinov and R. Kashapov, Theor. Math. Phys. 
183, 491 (2015). 

[15] A. G. Marinopoulos, L. Reining, A. Rubio, and V. Ole- 
vano, Phys. Rev. B 69, 245419 (2004). 

[16] J. F. Dobson, A. White, and A. Rubio, Phys. Rev. Lett. 
96, 073201 (2006); J. F. Dobson, K. McLennan, A. Ru¬ 
bio, J. Wang, T. Gould, H. M. Le, and B. P. Dinte, Aust. 
J. Ghem. 54, 513 (2001). 

[17] M. Bostrom and B. E. Sernelius, Phys. Rev. B 61, 2204 

( 2000 ). 

[18] M. G. Schabel and J. L. Martins, Phys. Rev. B 46, 7185 
(1992). 

[19] H. Rydberg, M. Dion, N. Jacobson, E. Schroder, 
P. Hyldgaard, S. 1. Simak, D. C. Langreth, and B. 1. 
Lundqvist, Phys. Rev. Lett. 91, 126402 (2003). 

[20] L. X. Benedict, N. G. Ghopra, M. L. Cohen, A. Zettl, 
S. G. Louie, and V. H. Crespi, Ghem. Phys. Lett. 286, 
490 (1998). 

[21] R. Zacharia, H. Ulbricht, and T. Hertel, Phys. Rev. B 
69, 155406 (2004). 


